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A thermodynamic model based on exergy flow through a Stirling Refrigerator is developed. Important 
irreversibilities of the refrigerator due to external heat transfer with the reservoirs, heat leak, flow and 
heat transfer in regenerator are included in the model. Expansion and compression efficiencies are intro¬ 
duced in the model to account for the losses in these processes. The effect of a control phase shift 
between the mass flow rate and pressure across regenerator on the performance of the refrigerator is pre¬ 
sented. Analytical solutions representing important quantities in the design of Stirling refrigerators such 
as the load curve, cooling power and efficiency in terms of basic system input parameters are developed. 
Thermodynamic bounds for the performance of Stirling refrigerators are obtained. Results indicating a 
compromise between cooling power and efficiency that are dependent on the constraint of the system 
are presented and discussed. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

Finite time thermodynamics and entropy generation minimiza¬ 
tion are established methods of performance analysis and optimi¬ 
zation of power and refrigeration cycles [1-3]. Different 
thermodynamic models of Stirling refrigerators, simplifying the 
complexity of cyclic heat and mass transfer in different compo¬ 
nents, have been subject of several studies [4-9], to just name a 
few. In this study a model for the exergy flow in Stirling refrigera¬ 
tor with emphasis on exergy analysis of regenerator in terms of its 
inlet and exit parameters is developed. In this manner, important 
sources of exergy destruction and dimensionless quantities con¬ 
trolling the cooling capacity and efficiency of the refrigerator are 
developed. Exergy analysis is a powerful method for the design 
of pulse tube and Stirling refrigerators [10-14]. Exergy flow and 
analysis of Stirling refrigerators show how the input exergy pro¬ 
vided by the power input to the compressor is destroyed as the 
working fluid goes through its cyclic motion in the system. In addi¬ 
tion, cooling capacity and efficiency are obtained in terms of 
important dimensionless numbers convenient for quantifying the 
important losses and performance evaluation of Stirling refrigera¬ 
tor at the system level. 

An important goal of this study is to find analytical expressions 
for the thermodynamic bounds of Stirling refrigerators in terms of 
the most important global variables influencing the performance of 
the refrigerator. Therefore, the study can be used to assess the ef- 
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feet that changing a global design parameter of the refrigerator can 
have on a particular aspect of its performance. Depending on the 
constraint on the system, the study shows that the efficiency and 
cooling power of the Stirling refrigerator can be two distinct crite¬ 
ria for the optimization of the refrigerator. Regenerator plays a par¬ 
ticularly important role on the thermodynamic bounds of Stirling 
refrigerators and greatly influences its cooling capacity and effi¬ 
ciency. For the purpose of this study, regenerator is conveniently 
modeled using two free parameters one representing its exergy 
flow efficiency and the other its thermal ineffectiveness. 

Exergy is a function of both the state of a system and the state of 
the environment as it measures the departure of the state of the 
system from the state of the environment. In application to Stirling 
refrigerators, for each component, considering one channel of heat 
transfer between the system and a thermal reservoir at the tem¬ 
perature T r and one channel of inlet and exit mass transfer, the 
exergy balance can be written as [10,12] 


Ed 


(lWe) r (Me) t - W + 



( 1 ) 


where e is the specific exergy carried with mass, E D is the rate of 
exergy destruction in the component, and Q is the rate of heat 
transfer. For a single-component working fluid, such as helium, in 
the absence of chemical, kinetic, or potential exergy, the specific 
exergy can simply be written as 

e = h-h 0 - T 0 (s - s 0 ) (2) 


where h is the specific enthalpy and s is the specific entropy and 
subscript “o” denotes the environment. Substituting Eq. (2) into 
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Nomenclature 



A 

area of regenerator 

X 

regenerator ineffectiveness 

C 

flow conductance 

CO 

angular frequency 

CHX 

Cold Heat Exchanger 

y 

specific heat ratio 

Cp 

heat capacity at constant pressure 

T 

period of oscillation 

Ed 

rate of exergy destruction 



e 

specific exergy 

Subscripts 

h 

specific enthalpy 

a 

average, ambient 

HHX 

Hot Heat Exchanger 

b 

bound 

/ 

rate of irreversibility 

c 

cold 

I< 

effective thermal conductivity of regenerator 

CO 

no load 

L 

length of regenerator 

cond 

conduction 

171 

amplitude of mass flow rate 

comp 

compressor 

M 

mass flow rate 

cn 

normalized cold side 

Mr 

mass flow rate amplitude ratio across regenerator 

cnl, cn2 

normalized cold side, Eqs. (24) and (29) 

P 

pressure 

e 

exit, expansion 

P 

pressure amplitude 

h 

hot heat exchanger 

Pr 

pressure amplitude ratio across regenerator 

ex 

exergetic 

Q 

rate of heat transfer 

i 

inlet 

s 

specific entropy 

i 

general inlet or exit of components 

T 

temperature 

i 

cold heat exchanger 

UA 

thermal conductance 

leak 

heat leak 

V 

average effective void volume of regenerator 

n 

normalized 

w 

power 

0 

environment 

Y 

sine of phase angles 

opt 

optimum 



p 

pressure 

Creeks 


r 

regenerator 

9 

phase angle for pressure 

R 

reservoir 

>1 

efficiency 

th 

thermal 


dimensionless number, m 2 /Cp } 

1 

hot side of regenerator, Fig. 1 

<1*2 

dimensionless number, Vco/C 

2 

cold side of regenerator, Fig. 1 

<t> 

phase angle for mass flow rate 

0, 1, 2, 3 

4 used with sine of phase angles, Appendix A. 


Eq. (1) and using the vanishing average mass flow rate over a cycle 
yields 

Ed = (M/i), - ( Mh) e - T 0 [(Ms), - (Ms)J -lV+(l-g<2 (3) 

Assuming the ideal gas law is valid and thermophysical proper¬ 
ties are constant, the specific enthalpy and entropy in the above 
equation can be written respectively as 

h = C„(T-T 0 ) (4) 

s= C p ln(T/r 0 ) - ftln(P/P 0 ) (5) 


2. Exergy based thermodynamic analysis of Stirling 
refrigerators 

Fig. 1 shows important components of a Stirling refrigerator 
used for analysis in this study. Exergy comes into the system at 
the driver side from the compressor and is destroyed as it moves 
through the system. The product exergy is delivered to the cold 
reservoir at the temperature T c . The conservation of energy and 
the balance of exergy for a system composed of compressor and 
the Hot Heat Exchanger (HHX) at the temperature T h can be writ¬ 
ten respectively as 



Fig. 1 . Schematic of Stirling refrigerator and parameters used in the model. 
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Qh — Wcomp + W e ~ Q r (6) 

+ VcompWcomp = E, + W = £, + Q* (l - Q (7) 

where W comp and W e are the input compressor power and power 
transferred from the expander over a cycle, respectively. £i is the 
rate of exergy transfer exiting the system, and Qh and Q r are the 
rate of heat transfer to the hot reservoir and energy transfer to 
regenerator, respectively. i HHX represent the rate of exergy destruc¬ 
tion in the HHX and f/ comp denotes the exergetic efficiency of com¬ 
pressor. We assume for purposes of this study that the mass flow 
rate and pressure at the inlet and exit of each component are given 
by 

Mj = rhj cos (cot + <f>j) (8) 

Pj = Pa+ Pj cos(cot + Bj) (9) 


Using Eqs. (2), (4), and (5), the pressure component of exergy at 
any location over a cycle can be written as [12] 



where t is the period of oscillation. Using Eqs. (8)—(10) the pressure 
component of exergy assuming relatively small pressure amplitude 
(Pj < 0.3 P„) can be written as 

(Ej) P = (l/2)J?T 0 rhj^ cos(0j - Bj) (11) 

1 a 

The thermal exergy transfer at the cold side of the regenerator 
due to conduction heat transfer and ineffectiveness of the regener¬ 
ator can be written as 


(£ 2 )* = Qt(i-^) (12) 

where T ( is the exit temperature of the regenerator. Using the defi¬ 
nition of thermal ineffectiveness of the regenerator, the rate of en¬ 
ergy transfer due to the mass flow at the cold side of regenerator 
can be estimated [12,15] 


(\/n)C p (T h m~T,m 2 ) 

where the denominator represents the maximum rate of enthalpy 
transfer into the regenerator. Conduction heat transfer to the cold 
side of regenerator can be estimated assuming a linear temperature 
profile in the regenerator 

Qcond = ( KA/L)(T h - T,) (14) 

where KA/L is effective thermal conductance in the regenerator and 
its shell. It should be pointed out that the direction of the thermal 
exergy transfer at the cold side of the regenerator is opposite to 
the direction of heat transfer at the point. Using Eqs. (13) and 
(14), the rate of heat transfer at the cold side of regenerator can 
be written as 


Qr = 0.2 + Qcond = (1 /n)XC p (T„m - T,ifi 2 ) + (I(A/L)(T h - T,) (15) 

Using Eqs. (11) and (12) the total exergy entering and leaving 
the regenerator can be written respectively as 

£i = (£i) P + (£i) tt 

= (1 /2)RT„m COS(^ - Bj) + Q r (l - Q (16) 

E 2 = (E 2 ) p + (E 2 ) th 

= (l/2)i?T 0 m 2 g cos(<)» 2 - fc) + Q r (l (17) 


Exergy balance for the Cold Heat Exchanger (CHX) can be writ¬ 
ten as 







(18) 


where rj e is the exergetic efficiency of the expander, Qieak is heat leak 
to CHX, and Q c and T c are the cooling capacity of the refrigerator 
and the cold reservoir temperature, respectively. The first and last 
terms of Eq. (18) represent the irreversibility of CHX and the exergy 
delivered to the cold reservoir, respectively. Using conservation of 
energy for CHX and expander, Eq. (18) can be used to find the cool¬ 
ing capacity of the refrigerator. 


Qc = 


'7c(£2) p 


> 1 e T 0 /T, -r] e + 1 


Qr Qieak 


(19) 


The heat transfer analysis of HHX, CHX and the heat leak is sim¬ 
plified using overall heat conductance for the heat exchangers and 
the heat leak. 


Q h = UA h (T h - T 0 ) (20) 

Qc = UAc(T c - Ti) (21) 


Qieak = UAi eak (T a - T t ) 


( 22 ) 


where UA is the effective thermal conductance and subscripts It, c 
and leak denote hot, cold and leak, respectively. It is assumed that 
the heat leak to CHX occurs from a reservoir at the effective temper¬ 
ature of T a . 

It is convenient to define the exergetic efficiency of flow through 
the regenerator as the ratio of the pressure component of exergy at 
its cold side to its hot side. Using the definition of the exergetic effi¬ 
ciency of flow through the regenerator Eq. (19) can be written as 


Qc = W„ 


>1e>lATo/Th + ri comp -V 
' 1 + HATo/T, - ipTo/Th - 1 ) 


Qr Qieak 


(23) 


where tj r is the exergetic efficiency of flow through the regenerator 
and can be obtained from the application of Eq. (11) to the inlet and 
exit of the regenerator. The above equation gives the cooling capac¬ 
ity of Stirling refrigerator in terms of important system parameters 
and exergetic efficiencies of its major components. For an ideal case 
were the heat conductance of the cold and hot heat exchangers are 
infinity, there is no heat leak, the regenerators thermal losses are 
zero, the efficiencies of the compressor, the expander, and the 
regenerator are one, Eq. (23) gives the Coefficient of Performance 
(COP), COP = Qc/Wcomp = T C /(T 0 - T c ), equal to the Carnot COP. In 
application of Eq. (23) the relations between important parameters 
of the system and the quantities W CO m P , Qr and t; r must be obtained. 

Substituting Eqs. (14), (15), and (22) into Eq. (23) and dividing 
the results by the quantity m^RT 0 yields 

A W compp PePrC^o/T h + ^7 comp 1) ^ y ( Th Mr Tl \ 

^ cnl ^+ri e (T 0 /T l -r lr T 0 /T h -\) ny~\\T 0 T 0 ) 

(■ KA/L) n (T„ - Ti) UA leakn (T a - T,) 

To T 0 

where Q cnl = Q c /miRT 0 , VV comPi „ = W comp /m,kT 0 , (I<A/L) n = (KA/L)/ 
theR, UAieahn = UAieait/rhjR, and Mr=m 2 /rhi. Therefore, Eq. (24) 
gives the normalized cooling capacity of the Stirling refrigerator 
in terms of important global system parameters. 

For practical application of the thermodynamic analysis pre¬ 
sented in this study, relations between the mass flow, pressure 
and phase angles at the cold and hot sides of regenerator must 
be obtained. Oscillating flow of heat and mass transfer in regener¬ 
ators is very complicated [15]. For thermodynamic analysis we 
consider a simple model for the regenerator to find the bounds 
for cooling capacity and efficiency of Stirling refrigerators. The 
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model considers the most important parameters affecting the 
regenerator losses including the effect of void in the regenerator. 
The first order model of the regenerator assumes that the mass 
flow rate at the cold side of the regenerator can be obtained using 
a linear relation based on a given appropriately average flow 
conductance. 

m 2 cos(cot + cj> 2 ) = C[p, cos(cot + 0i) - p 2 cos(cot+ 0 2 )] (25) 

where C is properly averaged regenerator flow conductance. In 
addition, a simple model for the conservation of mass in the regen¬ 
erator can be written as 

rh, cos(mt + <J> : ) = m 2 cos (cot + 4> 2 ) - Voj[p i sin(ftjt+ 0j) 

+ p 2 sin(cot + 0 2 )] (26) 

where V is the properly averaged effective regenerator void volume 
including the influence of temperature distribution in the regener¬ 
ator. Expanding Eqs. (25) and (26), and applying trigonometric 
identities, results in a set of equations that uniquely determine 
the relations between the phase shifts and amplitudes of pressure 
and mass flow rates at the cold and hot sides of the regenerator. 
These relations are given in Appendix A. Using Eq. (11 ), the exerget- 
ic flow efficiency of regenerator d r can be written as 


sume that the conduction heat transfer in the regenerator is zero, 
there is no heat leak, and no void in the regenerator, Eqs. (17), 
(19) and Al can be combined to find the cooling capacity of the 
refrigerator 


Qcn2 = 


>; c Prcos (<j) 2 - 0 2 ) _ 27 y Pa / TA' 
riJ 0 /T c -ri e +1 7ry-lp,V T 0 ) 

1 + Pr 2 - 2Prf \J cos 2 (<^> 2 - 0 2 )(1 - Pr 2 sin 2 (0 2 


fc)) 


+Prsin 2 (<£ 2 - 0 2 ))] 


1/2 


(29) 


where Q cn2 = 2P 0 Q/CRT„p^. Eq. (29) shows the normalized cooling 
capacity as a function of the pressure amplitude ratio and the phase 
shift between the mass flow and pressure at the cold side of the 
regenerator. The thermodynamic bound of the cooling capacity oc¬ 
curs when the phase shift is zero. In this case Eq. (29) can be consid¬ 
erably simplified and written as 


CLn2,b = (1 - Pr) 


lie p r 

JlJo/Tc - de + 1 


27 y Pa ( 1 _ Pc\ 

71 7 - 1 Pi V T 0 J_ 


(30) 


In this particular case the pressure amplitude ratio across the 
regenerator for the maximum cooling capacity can be analytically 
evaluated. 


;; r = Mr Pr 


cos(<fe - 0 2 ) 
cos(4>, - 0i) 


= Mr Pr 


1 -Yi 


1 -Yi 


(27) 


where Pr is the pressure amplitude ratio across the regenerator. 


3. Results and discussions 


Eq. (24) can be used to calculate the no-load temperature where 
the left hand side of the equation is zero. Using Eq. (21) this tem¬ 
perature is denoted by T co where T; = T c = T co . The no-load tem¬ 
perature is of great practical importance in design of cryogenic 
refrigerators [11], For the special case of infinite thermal conduc¬ 
tance of the HHX, no conduction heat transfer in the regenerator, 
no heat leak, and no void in the regenerator, Eq. (24) can be consid¬ 
erably simplified and the no-load temperature can be written as 
solution of a quadratic equation for T 0 /T C0 . 


i y T co ' 

Icomp'lAy-^K To, 



- w 


comp,n 


= o 


(28) 


Pf opt = 1/2 + A ^ (i _ Q(T„/T c + l/f7e - 1) (31) 

For the case of thermally perfect regenerator, when Pr opt = 1 /2, 
the thermodynamic bound for the normalized cooling capacity is 
given by Q cn2 , b = if 11 is interestin g t0 note that from 

Eqs. (30) and (31) for the Stirling refrigerator to be able to produce 
any cooling the upper bound of the regenerator ineffectiveness is 
given by 

7t(y-l) pi_1_ 

2y Pa (i-r c /r 0 )(T 0 /r c + i/f 7 e -i) 1 J 

Eqs. (30) and (31) can be used to find the thermodynamic 
bound for the no-load temperature for a given 7, p,/Pa and i) s . 

(1 - T C0 /T 0 )(t) e T 0 /T C0 -n e + 1) - nl7 2y lp a eP> = 0 (33) 

Assuming infinite thermal conductance for the hot and cold 
heat exchangers and using Eqs. (19) and (23), the expression for 
the compressor power can be written as 


where y is the specific heat ratio of the working fluid. For a given 
ratio of pressure amplitude to the system pressure at the inlet of 
the regenerator, the normalized compressor power is a slowly vary¬ 
ing function of the compressor power. Therefore, Eq. (28) shows the 
secondary effect of the compressor power and the significant influ¬ 
ence of the regenerator ineffectiveness on the no-load temperature 
of Stirling refrigerators. In addition, the equation also shows the ef¬ 
fect of exergetic efficiencies of the expander, compressor, and fluid 
flow through the regenerator on the no-load temperature. It should 
be pointed out that the reduction in the value of the regenerator 
ineffectiveness is normally accompanied by a reduction in the exer¬ 
getic efficiency of the fluid flow in the regenerator. This fact is well- 
known as an important challenge in the design of regenerators for 
the Stirling cryogenic refrigerators. 

3.1. Analysis for no void in regenerator with zero external 
irreversibility 

Thermodynamic analysis of Stirling refrigerators is considerably 
simplified when the thermal conductance of the cold and hot heat 
exchangers are infinity, resulting in zero exergy destruction in the 
heat exchangers (no external irreversibilities). In addition, if we as- 


_ [1 +^ e (T 0 /r c - ;/ r - 1 )]£ 2p 

'Ir'lcomp^eTo/Tc ~ de + 1) 


(34) 


Assuming no heat leak, Eqs. (23) and (34) can be combined to 
obtain an expression for the Coefficient of Performance (COP) of 
Stirling refrigerators 


COP = 


Q c 


dedrda 


Qr 


W a 


i + n e (T 0 /T c -f r -\) 


w„ 


(35) 


Using Eqs. (15) and (34) and assuming no conduction heat 
transfer in the regenerator, Eq. (35) can be written as 


COP = 


d ed compd r 


3 + >1e(To/T c — J/ r - 1) 

T-P ad pi comply (fl eT o ! T c ~ >1 e + 1)0 / Mr - T c /T 0 ) 

Tt Pi (y - 1)[1 + djTJTc - dr - l)]Prcos(</> 2 - 0 2 ) 


(36) 


In Eq. (36) )/ r is obtained from Eq. (27) and for the case of no 
void in the regenerator (Mr = 1) can be written as 


dr 


Prcos (4> 2 - 0 2 ) 

\J 1 - Pr 2 sin 2 (0 2 - 0 2 ) 


(37) 
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Exergetic efficiency can be obtained by dividing the COP by the 
coefficient of performance of the Carnot cycle, 
r] ex = COP/COP car = COP(T 0 /T c - 1). For the case of no void in the 
regenerator, Eqs. (36) and (37) indicate that the thermodynamic 
bound for the COP occurs when the mass flow and pressure are 
in phase at the cold side of regenerator. In this case the optimum 
COP can be written as 


rnp _ _ ^e^lcomp^ _ 

opt 1 + f] e (Jo/T c - Pr - 1) 

2P a r,c 0m M r lJo/Tc - lie + 1)0 - Tc/To) 

71 Pi(7 - 1)[1 + Vei^o/Tc - Pr - 1)] 1 j 

The COP corresponding to the optimum cooling capacity can be 
obtained by substituting Eq. (31) into Eq. (38). For a thermally per¬ 
fect regenerator, the optimum COP is given by the first term of Eq. 
(38). Eq. (31) gives Pr opt = 1/2 for a thermally perfect regenerator. 
In this case the corresponding coefficient of performance at opti¬ 
mum cooling capacity is given by COP opt = 1/(2 T 0 /T c - 1), also re¬ 
ported by de Boer [8]. It is convenient to write Eq. (36) for a given 
pressure phase shift between the hot and the cold sides of the 
regenerator. Eqs. (36) and (37) can be combined to find the COP 
in terms of the pressure phase shift across the regenerator utilizing 
the following trigonometric equation for the case of no void in the 
regenerator. 


tan(</> 2 - 0 2 ) 


sin(0] - 0 2 ) 
cos(0i - 0 2 ) - Pr 


(39) 


Therefore, for the case of thermally perfect regenerator the COP 
in terms of the pressure phase shift across the regenerator can be 
written as 


cop = _ yie'lc ompPr[cos(0i - 0 2 ) - Pr] _ 

1 - Prcos(0i - 0 2 ) + r/ e |j/-(l - Prcos(0i - 0 2 )) + Pr 1 2 - lj 

(40) 

The above equation is a convenient analytical expression to find 
the thermodynamic bound for the COP for a given pressure phase 
shift across the regenerator and different values of expander and 
compressor efficiencies. Based on Eq. (40), the thermodynamic 
bound occurs when the pressure phase shift between the hot and 
cold side of the regenerator is zero. Therefore, the optimum COP, 
including the effects of the compressor and expander efficiencies, 
can be written as 


COP 


opt 


lelc om p Pr(l - Pr) 

1 - Pr + r] e ^(1 - Pr) + Pr 2 - lj 


(41) 


In addition, for the ideal case when the compression and expan¬ 
sion efficiencies are one, Eq. (41) gives the coefficient of perfor¬ 
mance COP 0 pt = Pr /(T 0 /T c - Pr), also noted by de Boer [8]. 

Fig. 2 shows the results of calculations using Eq. (29) for the 
normalized cooling capacity as a function of the pressure ampli¬ 
tude ratio across the regenerator. The results are given for three 
values of ineffectiveness and three values of the phase shift be¬ 
tween the mass flow and pressure at the cold side of the regener¬ 
ator. Other values of parameters are given in the figure. It can 
easily be seen from Eq. (29) that, with the assumptions used to de¬ 
velop the equation, the same cooling capacity is obtained when the 
phase of mass flow is lagging or leading the pressure phase by the 
same amount. In addition, the thermodynamic bound for the cool¬ 
ing capacity occurs when the phase shift is zero as discussed pre¬ 
viously. The figure shows that the value of the regenerator 
ineffectiveness has a significant effect on the cooling capacity. In 
fact, extrapolating the values of the optimum cooling capacity for 
the zero phase shift, the upper bound for the regenerator ineffec- 



Fig. 2. Normalized cooling capacity as a function of Pr for different values of X and 
Vo- 


tiveness capable of producing any cooling is X « 0.0155. This value 
is consistent with the upper bound for the regenerator ineffective¬ 
ness given analytically by Eq. (32) for the same parameters used in 
the figure. Fig. 3 shows the same results as in Fig. 2 for different 
values of the regenerator ineffectiveness and the expander effi¬ 
ciency. The results are given only when the mass flow and pressure 
are in phase at the cold side of the regenerator. The optimum val¬ 
ues of Pr obtained from the figure correspond to the values ob¬ 
tained analytically from Eq. (31). The results also show the 
significant effect of the regenerator ineffectiveness on the cooling 
capacity of the refrigerator while the effect of the expander effi¬ 
ciency is not very significant for the values used in this figure. 
Fig. 4 shows exergetic efficiency as a function of the pressure 
amplitude ratio for two values of the regenerator ineffectiveness, 
different values of the phase shift at the cold side of regenerator, 
and different values of the pressure phase shift across the regener¬ 
ator. Other values used in the calculation are reported on the fig¬ 
ure. The results are qualitatively different when the phase shift 
constraint between the mass flow and pressure at the cold side is 
compared with the pressure phase shift constraint across the 
regenerator. Exergetic efficiency is a monotonically increasing 
function of Pr when Yo is given while there is an optimum when 
the Y1 is the constraint. The calculations in Fig. 4 are based on 
Eqs. (36), (37), and (39) with Pr as the independent variable. As 
in Fig. 2, the thermodynamic bound for the exergetic efficiency oc¬ 
curs when the pressure and mass flow are in phase at the cold side 



Fig. 3. Normalized cooling capacity as a function of Pr for different values of X and 

ti e . 
















236 


A. Razani et al./Cryogenics 50 (20i0) 231-238 



0 0 2 04 0.6 0.8 1 

Pr 


Fig. 4. Exergetic efficiency as a function of Pr for different values of A and Y 0 . 


parameter. The significant influence of the regenerator ineffective¬ 
ness is clearly shown in the figure. Fig. 6 shows the cooling capac¬ 
ity and efficiency diagram for different values of the pressure phase 
shift across the regenerator using Pr as a parameter. Calculations 
are based on Eqs. (29), (36), (37), and (39). In these calculations 
the values of V, are fixed and Pr is used as a parameter. The 
loop-shaped curves in the figure show the compromise between 
the cooling capacity and efficiency. The values of Pr for the opti¬ 
mum cooling capacity and optimum efficiency are different. As 
the value of the regenerator ineffectiveness increases the loop¬ 
shaped curve becomes smaller, indicating diminished performance 
of the refrigerator in both cooling capacity and efficiency. In fact, 
for the case of Y i = ±0.2 and X = 0.008 the performance of the 
refrigerator is drastically reduced for the parameters used in the 
figure. 


3.2. The effect of void in the regenerator and external irreversibilities 


of regenerator. In addition, the efficiency is the same when the 
mass flow is lagging or leading the pressure by the same amount. 
Fig. 5 shows the cooling capacity and efficiency diagram for differ¬ 
ent values of the regenerator ineffectiveness, expansion efficiency, 
and the phase shift between mass flow and pressure at the cold 
side of the regenerator. The figure is generated by using Pr as a 



Fig. 5. Normalized cooling capacity and efficiency diagram for different values of X 
and Vo. 



0 0.01 0.02 0.03 0.04 0.05 

Normalized cooling capacity 

Fig. 6. Normalized cooling capacity and efficiency diagram for different values of A 
and Y,. 


The void in the regenerator has important effect on the perfor¬ 
mance of the Stirling refrigerators. The dimensionless variable rep¬ 
resenting the void in the regenerator is given by i j/ 2 = Vco/C. The 
relations between the phase shifts at the inlet and exit of the 
regenerator, including the effect of void, is given in Appendix A. 
Fig. 7 shows the effect of the void in the regenerator on the normal¬ 
ized cooling capacity and the exergetic efficiency as a function of 
the pressure amplitude ratio across the regenerator for three val¬ 
ues of i jj 2 when the phase shift between the mass flow and pressure 
at the cold side of the regenerator is zero. The important effect of 
void in the regenerator is clearly shown. In addition, for one value 
of the void, \\i 2 = 0.2, the figure also gives the results for three val¬ 
ues of the phase shift between the mass flow and pressure at the 
cold side of regenerator. For a given void in the regenerator, the 
optimum cooling capacity and efficiency are different and are a 
function of the phase shift. Fig. 8 shows the cooling capacity and 
efficiency diagram for the same data used in Fig. 7 with Pr as a 
parameter. With Pr as a parameter, a looped-shape curve can be 
obtained indicating a compromise between the cooling capacity 
and efficiency for the cases studied. The optimum efficiency and 
the optimum cooling capacity are functions of the void in the 
regenerator and the phase shift as shown it the figure. For a given 
void in the regenerator, the optimum phase shift can be obtained 
for the cooling capacity or the efficiency. The figure indicates the 
potential advantage of actively controlling the phase shift in the 
Stirling refrigerator to optimize the cooling capacity or efficiency 
as compared to the passive methods of phase shift control in pulse 
tube refrigerators. 



Fig. 7. Normalized cooling capacity and exergetic efficiency as a function of the 
pressure amplitude ratio across regenerator for different values of i j/ 2 and Y 0 . 
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Fig. 8. Normalized cooling capacity and efficiency diagram for different values of 
i// 2 and V 0 . 

The effect of external irreversibilities due to the finite heat con¬ 
ductance of the heat exchange between the thermal reservoirs and 
the Stirling refrigerator is given in Fig. 9. In cryogenic refrigerators 
the magnitude of heat exchange in the hot heat exchanger is much 
larger than the cold heat exchanger. In this study the normalized 
thermal conductance of the heat exchangers is defined by 
UA n = 2(UA)P a /CRT 0 p Fig. 9 shows the cooling capacity and effi¬ 
ciency diagram for three equal values of the normalized thermal 
conductance for the hot and cold heat exchangers. As expected, 
the performance of the Stirling refrigerator is reduced when the 
irreversibility of heat exchange increases. In addition, the figure in¬ 
cludes the results for the case when the total heat conductance of 
the heat exchangers is fixed and the allocation of the thermal con¬ 
ductance to the heat exchangers is varied. The results show the 
well-known fact that an optimum exists for the allocation of the 
thermal conductance to the cold and the hot heat exchangers. 

3.3. Comparison with the Stirling cryocooler data 

In order to directly compare the experimental results for a par¬ 
ticular Stirling cryocooler to the control thermodynamic model 
used in this study, the global control parameters for the cryocooler 
must be known. Usually Stirling cryocoolers are designed for a 
specified set point and may not perform as well in off-design con¬ 
ditions; furthermore, global control parameters are not reported. 
The most important data usually reported for Stirling type cryoco- 



0 0.01 0.02 0.03 

Normalized cooling capacity 

Fig. 9. Normalized cooling capacity and efficiency diagram for different values of 
UA hn and UA m . 


olers are the load curves (Q c vs. T c ) for different ambient temper¬ 
atures and the compressor power. The no-load temperature and 
the coefficient of performance are also usually reported. As an 
example, the performance data for the CryoTel™ Stirling cryocool¬ 
er developed by Sunpower are selected to show that the results 
from the characterization of the cryocooler correspond to the re¬ 
sults obtained in this study [16,17]. One of the characteristics of 
the load curves for Stirling type cryocooler is the near linear vari¬ 
ation of load curves for a given compressor power. This can be eas¬ 
ily seen from Eq. (36) where for a given compressor power the 
variation of Q r is a nearly linear function of T c . In addition, for a gi¬ 
ven ambient temperature, the no-load temperature is a slowly 
varying function of the compressor power. This behavior can be 
seen from Eq. (28) where the combined effect of the thermody¬ 
namic parameters is a slowly varying function of the compressor 
power. For example, for the CryoTel CT Stirling cryocooler, a com¬ 
pressor power of 40 W produces a no-load temperature of approx¬ 
imately 52 I< while a compressor power of 160 W results in a no- 
load temperature of approximately 32 K [16]. This characteristic 
of the load curves is shown in the performance curves of many Stir¬ 
ling type cryocoolers characterized in our laboratory [11], The ef¬ 
fect of ambient temperature on the performance of Stirling 
cryocooler can be easily seen in analytical expressions developed 
in this study because To is a control parameter in the model and 
is also used as a parameter for characterization of Stirling crycool- 
ers. A recent single stage Sunpower Stirling cryocooler (CryoTel LT) 
can produce a very low no-load temperature of about 21.5 K for the 
ambient temperature of 310 I< [17]. It is interesting to compare this 
temperature to the thermodynamic bound developed in this study 
given by Eq. (33). As can be seen from the equation, this is only 
possible by designing a regenerator for the cryocooler with a very 
low value of thermal ineffectiveness. The regenerator for the Cryo¬ 
Tel LT single stage Stirling cryocooler is carefully designed to pro¬ 
duce such a low no-load temperature [17]. Eq. (36) clearly shows 
the important influence of compressor efficiency and the regener¬ 
ator exergetic flow efficiency along with its ineffectiveness on the 
efficiency of the Stirling refrigerator. Well-designed Stirling cryo¬ 
coolers currently reach the exergetic efficiency of 25% [16,17], Re¬ 
cently, the increase in compressor efficiency is expected to result in 
the development of a Stirling cryocooler with exergetic efficiency 
of close to 30% [18]. It should be pointed out that, as shown in this 
study, exergetic efficiency and cooling capacity can be two differ¬ 
ent optimization criteria in the design of Stirling crycoolers. 


4. Conclusions 

A thermodynamic model based on exergy flow for application 
to Stirling refrigerators is developed. The model includes the 
important sources of exergy destruction in the refrigerator. For 
the case of no void in the regenerator and no exergy destruction 
in the hot and cold heat exchangers, analytical solutions and ther¬ 
modynamic bounds for the performance of the refrigerator are ob¬ 
tained. The pressure ratio across the regenerator plays an 
important role in the development of the analytical solutions in 
this study. It is shown the under particular constraints one can ob¬ 
tain a loop-shaped curve for the cooling capacity and efficiency 
diagram, indicating a compromise between the cooling capacity 
and efficiency. The significant influence of the regenerator effec¬ 
tiveness and the void in the regenerator on the refrigerator perfor¬ 
mance is noted. An exergetic efficiency parameter defined for the 
fluid flow in the regenerator facilitates the thermodynamic analy¬ 
sis of the Stirling refrigerator and the development of thermody¬ 
namic bounds for the refrigerator. The effect of the size of hot 
and cold heat exchangers, represented by their values of thermal 
conductance on the performance of the Stirling refrigerator, is eval- 
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uated. In this study we have formulated the results in terms of 
non-dimensionalized variables which are convenient for optimiza¬ 
tion, parametric studies and obtaining the thermodynamic bounds 
of the Stirling refrigerators. 


eters i/q, Mr, Y, = sin(0i - 0 2 ), Y 2 = sin(^> 1 - 0 2 ). Y 3 = sin^, - rii), 
and Y 4 = sin(0, - <f> 2 ) can be obtained. 
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Appendix A 


Important phase angles for the pressure and mass flow rate 
across regenerator in terms of dimensionless parameters are given 
below. These equations are the results of the trigonometric identi¬ 
ties applied to Eqs. (25) and (26). Therefore, all trigonometric func¬ 
tions used in Eqs. (25) and (26) can be uniquely determined. 


l/q = 


1 + Pr 2 - 2Pr^ v /(l - Y 2 )(l -Pr 2 Y 2 ) + PrY 2 ^) 


1/2 


(Al) 


Yi =Yn \/l — Pr Y„ — Pr\/1 - Y; 


(A2) 


Mr 2 + 1 - 2Mr sj{\ - Y 2 )(l - Y 2 ) + Y 0 Y 2 
- (Mri/f 2 /i/q) 2 (l + Pr 2 + 2Pn/l - Y 2 j = 0 

s](\- Y 2 )(l - Y 2 ) + Y,Y 2 - Mr [7(1 - Y 2 )(l - Y 2 ) + Y 0 Y, 
+ Pru / 1 - Y 2 - MrVrJ 1 - Y 2 = 0 


(A3) 


(A4) 


Y 3 = VW1 -Y\~yJ\ -Y 2 


(A5) 


where the dimensionless parameters i/q = m 2 /Cp^ and \\i 2 = Vco/C 
are defined representing the physical and thermal properties of 
the regenerator. Therefore, using the above five equations, for se¬ 
lected parameters such as Y 0 = sin(0 2 - 0 2 ), Pr, and 1 l/ 2 , the param- 
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